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The method of the nonequilibrium statistical operator developed by D. N. Zubarev is employed to 
■ analyse and derive generalized transport and kinetic equations. The degrees of freedom in solids can 

D ' often be represented as a few interacting subsystems (electrons, spins, phonons, nuclear spins, etc.). 

P-h , Perturbation of one subsystem may produce a nonequilibrium state which is then relaxed to an 

equilibrium state due to the interaction between particles or with a thermal bath. The generalized 
kinetic equations were derived for a system weakly coupled to a thermal bath to elucidate the 
nature of transport and relaxation processes. It was shown that the "collision term" had the same 
functional form as for the generalized kinetic equations for the system with small interactions among 
O ' particles. The applicability of the general formalism to physically relevant situations is investigated. 

Sit is shown that some known generalized kinetic equations (e.g. kinetic equation for magnons, 
Peierls equation for phonons) naturally emerges within the NSO formalism. The relaxation of a 
small dynamic subsystem in contact with a thermal bath is considered on the basis of the derived 
equations. The Schrodinger-type equation for the average amplitude describing the energy shift and 
damping of a particle in a thermal bath and the coupled kinetic equation describing the dynamic 
and statistical aspects of the motion are derived and analysed. The equations derived can help in 
the understanding of the origin of irreversible behavior in quantum phenomena. 
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I. INTRODUCTION 



> 

3^ . The aim of statistical mechanics is to give a consistent formalism for a microscopic description of macroscopic 
behavior of matter in bulk. The methods of equilibrium and nonequilibrium statistical mechanics have been fruitfully 
, applied to a large variety of phenomena and materials Q, H H El H H • The statistical mechanics of irreversible 
f^i ■ processes in solids, liquids, and complex materials like a soft matter are at the present time of much interest. The 
lO | central problem of nonequilibrium statistical mechanics is to derive a set of equations which describe irreversible 
processes from the reversible equations of motion. The consistent calculation of transport coefficients is of particular 
interest because one can get information on the microscopic structure of the condensed matter. There exist a lot 
of theoretical methods for calculation of transport coefficients as a rule having a fairly restricted range of validity 
and applicability. The most extensively developed theory of transport processes is that based on the Boltzmann 
equation H |j| . However, this approach has strong restrictions and can reasonably be applied to a strongly rarefied 
gas of point particles (Tfl- For systems in the state of statistical equilibrium, there is the Gibbs distribution by 
means of which it is possible to calculate an average value of any dynamical quantity. No such universal distribution 
has been formulated for irreversible processes. Thus, to proceed to the solution of problems of statistical mechanics 
of nonequilibrium systems, it is necessary to resort to various approximate methods 0, flU. ILH ITEl fl6l ITtI Il8| . 
Kubo and others |19l 1 2(1 121| derived the quantum statistical expressions for transport coefficients such as electric and 
thermal conductivities. They considered the case of mechanical disturbances such as an electric field. The mechanical 
disturbance is expressed as a definite perturbing Hamiltonian and the deviation from equilibrium caused by it can be 
obtained by perturbation theory. On the other hand, thermal disturbances such as density and temperature gradients 
cannot be expressed as a perturbing Hamiltonian in an unambiguous way. During the last decades, a number of 
schemes have been concerned with a more general and consistent approach to transport theory 0, l22l I23L I24I I25L 
in mm in 113. These approaches, each in its own way, lead us to substantial advances in the understanding of 
the nonequilibrium behavior of many-particle classical and quantum systems. In addition, they have used dynamic 
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arguments to obtain kinetic and balance equations which describe the irreversible evolution of a system from particular 
initial states. This field is very active and there are many aspects to the problem |3l| . The purpose of the present 
work is to elucidate further the nature of transport processes and irreversible phenomena from a dynamic point of 
view. According to Montroll |32, "dynamics is the science of cleverly applying the operator exp(—iHt/h)". We wish 
to give a self-contained consideration of some general approach to the description of transport phenomena starting 
with dynamic equations. Our purpose here is to discuss the derivation, within the formalism of the nonequilibrium 
statistical operator p2l I27I l33l ] . of the generalized transport and kinetic equations. On this basis we shall derive, by 
statistical mechanics methods, the kinetic equations for a system weakly coupled to a thermal bath. 
In section [HI we briefly review some basic concepts. In section III Al the derivation of the transport and kinetic 
equations within the NSO formalism is outlined. In section 111 Jr?l we consider the application of the established 
equations to the derivation of the kinetic equations for magnons and phonons. Special attention is given to the 
problem of derivation of kinetic equations for a system weakly coupled to a thermal bath in section ITTT1 On the basis 
of these equations the balance and master equations are obtained in section llVl The behavior of a small dynamic 
system weakly coupled to a thermal bath is discussed in some detail in section The relaxation of a small dynamic 
subsystem in contact with a thermal bath is considered on the basis of the derived equations. The Schrodinger-type 
equation for an average amplitude describing the energy shift and damping of a particle in a thermal bath, and the 
coupled kinetic equation describing the dynamic and statistical aspects of the motion are derived and analysed in 
section IVll 



II. OUTLINE OF THE NONEQUILIBRIUM STATISTICAL OPERATOR METHOD 

In this section, we briefly recapitulate the main ideas of the nonequilibrium statistical operator approach [23 . 12H l33| 
for the sake of a self-contained formulation. The central statement of the statistical-mechanical picture is the fact 
that it is practically impossible to give a complete description of the state of a complex macroscopic system. We 
must substantially reduce the number of variables and confine ourselves to the description of the system which is 
considerably less then complete. The problem of predicting probable behavior of a system at some specified time is a 
statistical one. As it was shown by Gibbs an d Boltzmann it is useful and workable to employ the technique 
of representing the system by means of an ensemble consisting of a large number of identical copies of a single system 
under consideration. The state of the ensemble is then described by a distribution function p(fi . . . r n ,pi . . .p n , t) in 
the phase space of a single system. This distribution function is chosen so that averages over the ensemble are in exact 
agreement with the incomplete ( macroscopic ) knowledge of the state of the system at some specified time. Then 
the expected development of the system at subsequent times is modelled via the average behavior of members of the 
representative ensemble. It is evident that there are many different ways in which an ensemble could be constructed. 
As a result, the basic notion, the distribution function p is not uniquely defined. Moreover, contrary to the description 
of a system in the state of thermodynamic equilibrium which is only one for fixed values of volume, energy, particle 
number, etc., the number of nonequilibrium states is large. The role of the relaxation times to equilibrium state 
was analysed in paper (34[. The precise definition of the nonequilibrium state is quite difficult and complicated, and 
is not uniquely specified. Since it is virtually impossible and impractical to try to describe in detail the state of a 
complex macroscopic system in the nonequilibrium state, the method of reducing the number of relevant variables 
was widely used. A large and important class of transport processes can reasonably be modelled in terms of a reduced 
number of macroscopic relevant variables [3^ . There are different time scales and different sets of the relevant vari- 
ables [3(||37|], e.g. hydrodynamic, kinetic, etc. This line of reasoning has led to seminal ideas on the construction of 
Gibbs- type ensembles for nonequilibrium systems |28l l38l l39l l40j . B. Robertson |4lLl42l I43L l44j proposed the method 
of equations of motion for the "relevant" variables, the space- and time-dependent thermodynamic "coordinates" of a 
many-body nonequilibrium system which were derived directly from the Liouville equation. This was done by defining 
a generalized canonical density operator depending only upon present values of the thermodynamic " coordinates" . 
The most satisfactory and workable approach to the construction of Gibbs-type ensembles for the nonequilibrium 
systems, as it appears to the writer, is the method of nonequilibrium statistical operator (NSO) developed by D. N. 
Zubarev [27]] • The NSO method permits one to generalize the Gibbs ensemble method to the nonequilibrium 
case and to construct a nonequilibrium statistical operator which enables one to obtain the transport equations and 
calculate the kinetic coefficients in terms of correlation functions, and which, in the case of equilibrium, goes over to 
the Gibbs distribution. Although this method is well known, we shall briefly recall it, mostly in order to introduce 
the notation needed in the following. 

The NSO method sets out as follows. The irreversible processes which can be considered as a reaction of a system 
on mechanical perturbations can be analysed by means of the method of linear reaction on the external perturba- 
tion [l^. However, there is also a class of irreversible processes induced by thermal perturbations due to the internal 
inhomogeneity of a system. Among them we have, e.g., diffusion, thermal conductivity, and viscosity. In certain 
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approximate schemes it is possible to express such processes by mechanical perturbations which artificially induce 
similar nonequilibrium processes. However, the fact is that the division of perturbations into mechanical and thermal 
ones is reasonable in the linear approximation only. In the higher approximations in the perturbation, mechanical 
perturbations can lead effectively to the appearance of thermal perturbations. 

The NSO method permits one to formulate a workable scheme for description of the statistical mechanics of irreversible 
processes which include the thermal perturbation in a unified and coherent fashion. To perform this, it is necessary to 
construct statistical ensembles representing the macroscopic conditions determining the system. Such a formulation 
is quite reasonable if we consider our system for a suitable large time. For these large times the particular properties 
of the initial state of the system are irrelevant and the relevant number of variables necessary for description of the 
system reduces substantially [3(|. 

As an introduction to the NSO method, let us describe the main ideas of this approach as follows. The basic hypoth- 
esis is that after small time-interval r the nonequilibrium distribution is established. Moreover, it is supposed that it 
is weakly time-dependent by means of its parameter only. Then the statistical operator p for t > r can be considered 
as an "integral of motion" of the quantum Liouville equation 

Here |j| denotes time differentiation with respect to the time variable on which the relevant parameters F m depend. 
It is important to note once again that p depends on t by means of F m (t) only. We may consider that the system is in 
thermal, material, and mechanical contact with a combination of thermal baths and reservoirs maintaining the given 
distribution of parameters F m . For example, it can be the densities of energy, momentum, and particle number for the 
system which is macroscopically defined by given fields of temperature, chemical potential and velocity. It is assumed 
that the chosen set of parameters is sufficient to characterize macroscopically the state of the system. The set of the 
relevant parameters are dictated by the external conditions for the system under consideration and, therefore, the 
term ^ appears as the result of the external influence upon the system. Due to this influence precisely, the behavior 
of the system is nonstationary. 

In order to describe the nonequilibrium process, it is necessary also to choose the reduced set of relevant operators 
P m , where m is the index ( continuous or discrete). In the quantum case, all operators are considered to be in the 
Heisenberg representation 

P m (t) = exp(^)P m exp(^^) (2) 
n n 

where H does not depend on the time. The relevant operators may be scalars or vectors. The equations of motions 
for P m will lead to the suitable "evolution equations" .27]. In the quantum case 

«^a + ±[P. (( , lfl] -o. (3> 

The time argument of the operator P m (t) denotes the Heisenberg representation with the Hamiltonian H independent 
of time. Then we suppose that the state of the ensemble is described by a nonequilibrium statistical operator which 
is a functional of P m (t) 



p(t)=p{...P m (t)...} (4) 

Then p(t) satisfies the Liouville equation (JTJ. Hence the quasi-equilibrium ( "local-equilibrium") Gibbs-type distri- 
bution will have the form 



Pq = Q- 1 exp (- Fm(t)P m J (5) 

where the parameters F m (t) have the sense of time-dependent thermodynamic parameters, e.g., of temperature, 
chemical potential, and velocity ( for the hydrodynamic stage), or the occupation numbers of one-particle states (for 
the kinetic stage). The statistical functional Q q is defined by demanding that the operator p q be normalized and 
equal to 



Q q = Tr exp ( - F ™ W P »> ) ( 6 ) 
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This description is still very simplified. There are various effects which can make the picture more complicated. The 
quasi-equilibrium distribution is not necessarily close to the stationary stable state. There exists another, completely 
independent method for choosing a suitable quasi-equilibrium distribution 0, El HE EH EH ■ For the state with the 
extremal value of the informational entropy 0, 

S = -Tr(p]np), (7) 

provided that 

Tr(pP m ) =< P m > q - Trp = 1, (8) 

it is possible to construct a suitable quasi-equilibrium ensemble. Then the corresponding quasi-equilibrium ( or local 
equilibrium ) distribution has the form 

Pq = exp I n - V F m (t)P m ) = exp{S(t, 0)) (9) 



exp m - F m (t)P m ^j = cxp(S(t, 0)) 
n = In Tr exp ^-^F m (t)P ro j 



where S(t, 0) can be called the entropy operator. The form of the quasi-equilibrium statistical operator was constructed 
in so as to ensure that the thermodynamic equalities for the relevant parameters F m (t) 

SIuQq ^ <P m >,, ^—-F m (t) (io) 



6F m (t) SF m (t) «' 5<P m > 



are satisfied. It is clear that the variables F m (t) and < P m > q are thermodynamically conjugate. Here the notation 
used is < . . . > q = Tr(p q . . .). 

It is clear, however, that the operator p q itself does not satisfy the Liouville equation. The quasi-equilibrium operator 
should be modified in such a way that the resulting statistical operator satisfies the Liouville equation. This is the 
most delicate and subtle point of the whole method. 

By definition a special set of operators should be constructed which depends on the time through the parameters 
F m (t) by taking the invariant part of the operators F m (t)P m occurring in the logarithm of the quasi-equilibrium 
distribution, i.e., 

fO 

B m {t) = F m (t)P m = s / e £tl F m (i + fi)P m (fi)dti = (11) 

J — OO 

F m (t)P m ~ f dhe eti (F m {t + ii)P m (ti) + F m (t + ii)P m (ii)) 



where (e — > 0) and 



in at 



The parameter e > will be set equal to zero, but only after the thermodynamic limit has been taken. Thus, the 
invariant part is taken with respect to the motion with Hamiltonian H. The operators B m (t) satisfy the Liouville 
equation in the limit (e — > 0) 



dRrr, 

dt 



i r° i 

— [B m ,H]=e / dhe £tl [F m (t + t 1 )P m (t 1 ) + F m (t + t 1 )P m (t 1 ) 

III J — OO ^ 



(12) 



The operation of taking the invariant part, of smoothing the oscillating terms, is used in the formal theory of scat- 
tering [47] | to set the boundary conditions which exclude the advanced solutions of the Schrodinger equation |48j . It 
is most clearly seen when the parameters F m (t) are independent of time. 
Differentiating P m with respect to time gives 

^^=e f e^P m {t + t x )dh (13) 
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The P m (t) will be called the integrals ( or quasi- integrals ) of motion, although they are conserved only in the limit 
(e — > 0). It is clear that for the Schrodinger equation such a procedure excludes the advanced solutions by choosing 
the initial conditions. In the present context this procedure leads to the selection of the retarded solutions of the 
Liouville equation. This philosophy has been pressed by the necessity of a consistent description of the irreversibility 
which is, according to [ijj, " at once a profound and an elusive concept" ( c.f., a discussion in Refs. |3ll l50|] 1 . 
It should be noted that the same calculations can also be made with a deeper concept, the methods of quasi- 
averages |27l l33i IKlj . Let us note once again that the quantum Liouville equation, like the classical one, is symmetric 
under time-reversal transformation. However, the solution of the Liouville equation is unstable with respect to 
small perturbations violating this symmetry of the equation. Indeed, let us consider the Liouville equation with an 
infinitesimally small source into the right-hand side 

?£l + ±-{ Pe ,H} = -e(p s -p q ) (14) 

or equivalently 

+ ^Pe,H] = -e(]np E -]np g ), (15) 

where (e — * 0) after the thermodynamic limit. This equation (|14J) is analogous to the corresponding equation of the 
quantum scattering theory |47l |48| . The introduction of infinitesimally small sources into the Liouville equation is 
equivalent to the boundary condition 

e (^i) {p(t + h) - p q (t + tl )) e (Z ^) - 0, (16) 
where t\ — oo after the thermodynamic limiting process. It was shown |27ll33j that the operator p E has the form 

p e (t,t)=e d*ie £ ( {l *i) = e / dt ie etl p q (t + t u t + h) (17) 

J — oo J — oo 

Here the first argument of p(t, t) is due to the indirect time-dependence via the parameters F m (t) and the second one 
is due to the Heisenberg representation. The required nonequilibrium statistical operator is defined as 

r° 

p e = p e (t,0)=p q (t,0)=e dhe^pgit + tuh) (18) 

J — oo 

Hence the nonequilibrium statistical operator can then be written in the form 



Q- 1 exp I - F m (t)P m +Y.J 



p = Q- 1 exp | ->> m ) =Q- i cxp I -}_^e I dhe stl (F m (t + *i)P m (ti)) = (19) 

dhe £tl [F m (t + h)P m (ti) + F m (t + ti)P m (*i)] ) 



Let us write down Eq. l|15|l in the following form: 

j t (e 6t \np(t,t)) =ee £t lnp 9 (M), (20) 

where 

jjjf 

In p(t,t) =UUt,0)lnp(t > 0)U(t,0); U(t, 0) = expf— ) (21) 

n 

After integration, Ea. (|20|l . over the interval (— oo,0) we get 

r° 

lnp(t,t)=e dt x e stl In p q {t + t u t + h) (22) 

J — oo 

Here we suppose that lim £ _ >0 + In p(t, t) = 0. 

Now we can rewrite the nonequilibrium statistical operator in the following useful form: 



p{t, 0) = exp (-e [ 
V J— i 



dhe^hipqit + h,^)) = exp(lnp g (f,0)) = exp 0)) (23) 
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The average value of any dynamic variable A is given by 

<A>= lim Tr(p(t,0)A) (24) 

and is, in fact, the quasi-average. The normalization of the quasi-equilibrium distribution p q will persists after taking 
the invariant part if the following conditions are required 

Tr{p{t, 0)P m ) =< P m >=< P m > q \ Trp = 1 (25) 

Before closing this section, we shall mention some modification of the "canonical" NSO method which was proposed 
in and which one has to take into account in a more accurate treatment of transport processes. 



A. The Transport and Kinetic Equations 

It is well known that the kinetic equations are of great interest in the theory of transport processes. Indeed, as 
it was shown in the preceding section, the main quantities involved are the following thermodynamically conjugate 
values: 

oF m (t) 6 < P m > 

The generalized transport equations which describe the time evolution of variables < P m > and F m follow from the 
equation of motion for the P m , averaged with the nonequilibrium statistical operator (|23() . It reads 

The entropy production has the form 

S(t) =< S(t, 0)>=-J2< p m > F m {t) = - £ AF f ? F F n (t)F m (t) (28) 

These equations are the mutually conjugate and with Eq.J2BJ> form a complete system of equations for the calculation 
of values < P rn > and F m . 

Let us illustrate the NSO method by considering the derivation of kinetic equations for a system of weakly interacting 
particles In this case the Hamiltonian can be written in the form 

H = H Q + V, (29) 

where Hq is the Hamiltonian of noninteracting particles ( or quasiparticles ) and V is the operator describing the 
weak interaction among them. Let us choose the set of operators P m — Pk whose average values correspond to 
the particle distribution functions, e.g., a^ak or a^ak+q- Here at and are the creation and annihilation second 
quantized operators ( Bose or Fermi type). These operators obey the following quantum equation of motion: 

Pk = ^[P k ,H] (30) 
in 

It is reasonable to assume that the following relation is fulfilled 

[Pk,H ] =Y, c klPu (31) 

l 

where c^i are some coefficients ( c- numbers). 

According to Eci. (|19|l . the nonequilibrium statistical operator has the form 



= Q- 1 expl-Y / Fk(t)Pk+Y / [ dhe^lFkit + t^Pkit^+Fkit + t^Pkih)} 



(32) 
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After elimination of the time-derivatives with the help of the equation < P k >=< P k > q it can be shown [53 that 
the integral term in the exponent, Eq. l|32[) . will be proportional to the interaction V. The averaging of Ea. H3U|l with 
NSO 1)32(1 gives the generalized kinetic equations for < > 

d < P k > 1 _ , 1 v-^ 1 . , 

= — < [P k , H] >= — V c kl < P l> +- < [P k , V] > 33 

in iti — ^ in 



dt ih ' ifi 



Hence the calculation of the r.h.s. of l(33|l leads to the explicit expressions for the "collision integral" ( collision terms). 
Since the interaction is small, it is possible to rewrite Eq. l|33|) in the following form: 



d< P k > 



dt 



Ll + Ll + Lf+Lf, (34) 



where 



i% 
i 



L\ = ^<[P k ,V]> q (36) 



L f = ^l dhe^ <[V(t 1 ),[P k ,V]]> q 
n J-oo 

Lf = ±f° dt 1 e^<[V(t 1 ), i h^P l dLl{ --- <Pl> --- )] 

<l J — GO i 



d<P l > 

The higher order terms proportional to the V 3 , V A , etc., can be derived straightforwardly. 



(37) 
(38) 



B. Kinetic Equations for Magnons and Phonons 

The dynamic behavior of charge > magnetic [U > an d lattice |H3 systems is of interest for the study of transport 
processes in solids. Partial emphasis has been placed on the derivation of the kinetic equations describing the hot 
electron transport in semiconductors [56l l57j , and the relaxation of magnons [Eli. l59| and phonons [55l l60l | due to the 
inelastic scattering of quasiparticles. 

We discuss briefly in this section the processes occurring after the switching off the external magnetic field in a ferro- 
magnetic crystal. Our main interest is in ferromagnetic insulators, where the dominant interaction is the Heisenberg 
exchange coupling JSiSj. It is well known that a strong microwave magnetic field applied parallel to the dc field can 
give rise to parametric excitation of spin waves [5^, 0, • In this technique the wave number of the potentially 
unstable spin waves can be changed by varying the dc magnetic field. One thus obtains information about the varia- 
tion of the spin- wave relaxation time with the wave number |63l |6J, l65l l66l l67| . Because of its relative simplicity the 
"parallel pumping" technique has proved very useful in determining rather fundamental properties of ferromagnetic 
materials. The subharmonic generation of spin waves at high power levels is an efficient research tool for probing 
magnon-magnon and magnon-phonon interactions. Useful information about the spin-wave relaxation rate can be 
deduced from the kinetic equations to study magnon-magnon and magnon-phonon interactions. 

Here the spin-wave relaxation processes arising from the dipolar interaction will be considered as an example. The 
Hamiltonian has the form 

H = — J] J{Rw)SiS v +2(M)h J2 S i + 2 t4J2^JT (RwSiSi' -3(IUi>Si>)(RwSi')) (39) 

This Hamiltonian contains Zeeman energy, exchange energy, and dipolar en erg y. To treat this Hamiltonian, it should 
be expressed in terms of the amplitudes of the normal modes or spin waves [681 ] . The amplitudes of the normal modes 
are quantum-mechanically interpreted as creation and annihilation operators ( usually bosons). We get |6S| 



S+ = Sf+iS? = V2Sb\\ll- b -^; Sf = Sf-iSy = V2S\ll-°-£b l ; s f = -S + hjI„ (-MM 
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We adopt the notation 

bi = N~ 1/2 h exp(ikRi), b\ = iV~ 1/2 Y h l exp(-ikRi) 

k k 

The transformed Hamiltonian contains a term that is quadratic in the spin-wave amplitudes H ^ and also terms that 
are of higher order, H^ 3 \ etc. 

H = i? (2) + H {3) + H {A) + ... (41) 

The eigenstates of the quadratic part of the Hamiltonian can be characterized by the occupation numbers c k c k , 
i.e., the quadratic part can be diagonalized to the form |68j ] 

= Y <k)c{c k - e(k) = huj(k) (42) 

k 

Here the operators ct and c k are the second-quantized operators of creation and annihilation of magnons. All higher 
order terms in the Hamiltonian lead to transitions between the eigenstates. In terms of the magnon operators these 
terms are given by 

=YHk,p,p')4clc p ,A(k + p-p>)+H.C. (43) 

kpp' 

ff (4) = Y ®(k,P;p',r)clclc p ,c r A(k+p-p> -r)+H.C. (44) 

kpp'r 

Usually, only that term in the Hamiltonian which is of the third order in the amplitudes of the normal modes is 
considered explicitly, because only this term leads to the relaxation rates proportional to the temperature in the 
high-temperature limit. 

Let us apply now the formalism of generalized kinetic equations, as described above. We suppose that the set of 
averages < P k >=< c\ck >~< n k > characterize the nonequilibrium state of the system. The quasi-equilibrium 
statistical operator has the form 



Pq = Q- 1 exp ^- Y Fk(t)n k ^J ; Q' 1 = Tr exp ^- £ F k {t)n k ^j (45) 



The kinetic equation l|34[) can then be expressed by 

d < n k > 



= Ll + Ll + Lf+Lf (46) 



dt 

For both the contributions and the following equality holds: 

L° k = L\ = Lf = (47) 
Let us first consider the contribution of the term . We can then write 

Lf = —YA^'Prf)? 6 ^ + ^ - +p -p>) (48) 

p,p' 

[(< n k > +1)(< n p > +1) < n p , > - < n k >< n p > (< ry > +1)] 
-~\Hk,p,p')\ 2 5(u>(k)-w(p)-u;(p'))A(k-p-P) 
[(< n k > +1) < n p >< n p , > - < n k > (< n p > +1)(< n p , > +1)] 

We can make the same calculation to obtain L 2 k for the magnon-magnon scattering term 

[(< n k > +!)(< n p > +1) < n p , >< n r > - < n k >< n p > (< ry > +!)(< n r > +1)] 
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Here the notation was used 

< n k >= N{huj(k)) = [exp(/3huj(k)) - I]" 1 

The quantities &(k,p,p') and <&{k,pjJ,r) are the combination of the matrix elements which describe the various 
transitions between spin eigenstates |58j . Equation 148|) corresponds precisely to the rate equation which describes 
the change of the average occupation number < n k > of the mode k derived in [5^|. The discussion of the two 
relevant relaxation rates tZ ( due to the confluence and splitting ) is given there. The types of kinetic equations, 
Eas. l48|l . 149|l . involved in our derivation and the conclusions arrived at show very clearly that the NSO method is 
a workable and useful approach for derivation of the kinetic equations for concrete physical problems. As far as the 
kinetic equations for magnons is concerned, its convenience can become even more evident if one needs to take into 
account higher order magnon processes ( four, five, etc.). The higher order processes may give rise to additional and 
unusual behavior ( i.e., a general heating of the spin- wave system causing the saturation, additional smaller peaks 
and kinks in the measured curves, etc.) 

It is evident that a similar derivation can be given for the kinetic equation for phonons. The theory of thermal 
conductivity E3 nas been extensively developed beginnin g w ith the kinetic theory of Peierls 0, IHfil l70l | . The 
theory of lattice thermal conductivity invented by Peierls [HEl l7(1| is based on the assumption that the perturbing 
mechanisms to the harmonic case ( anharmonicity, imperfections) are small in magnitude. The Peierls collision term 
for the three-phonon processes looks like 

L f ~ Y ^2{mk,p,p')\ 2 6 (cj(k) + W (p) - W (p0) (50) 
p,p' 

[(< n k > +1)(< n. p > +1) < ry > - < n k >< n p > (< ry > +1)] 
+^\<f>(k,p,p)\ 2 5 (w(fc) - lu( P ) - uj{p)) [(< n k > +1) < n p >< ry >] 

Note that our calculations show that three- and four-phonon processes behave quite differently. One expects that the 
stronger the anharmonicity the larger the thermal resistance. To catch this trend, some sophisticated formalisms [TH ] 
have been developed which utilize a modified version of the Peierls-Boltzmann equation. A useful approach to 
improving the initial Peierls theory corresponds to the derivation of a generalized Peierls-Boltzmann equation, where 
the phonons in the collision term are treated not as free phonons but as quasiparticles with a finite width and 
damping which are determined sclf-consistently. Crystal lattices at low temperatures represent an interacting system 
of quasiparticles in which we observe two relaxation mechanisms of widely different time scales, i.e., the system either 
at short or long times after its initial perturbation from equilibrium. For the long-time behavior of the system it is 
possible to formulate the problem in terms of the correlation functions of quantities relaxing slowly, such as densities 
of conserved variables in the system. The corresponding transport equations are similar in structure to the phonon 
Boltzmann equation with a modification of the collision term. A detailed study of the transport equations for phonon 
systems is not within the scope of this paper and deserves a separate consideration. 



III. SYSTEM IN THERMAL BATH: GENERALIZED KINETIC EQUATIONS 

We now proceed to derive generalized kinetic equations for the system weakly coupled to a thermal bath. Examples 
of such systems can be an atomic ( or molecular) system interacting with the electromagnetic field it generates as 
with a thermal bath, a system of electrons or exitons interacting with the phonon field, etc. Our aim is to investigate 
relaxation processes in two weakly interacting subsystems, one of which is in the nonequilibrium state and the other 
is considered as a thermal bath. The concept of thermal bath or heat reservoir, i.e., a system that has effectively 
an infinite number of degrees of freedom, was not formulated precisely. A standard definition of the thermal bath is 
a heat reservoir defining a temperature of the system environment. From a mathematical point of view [3^ |. a heat 
bath is somethin g th at gives a stochastic influence on the system under consideration. In this sense, the generalized 
master equation |72l 1731 i s a tool for extracting the dynamics of a subsystem of a larger system by the use of a special 
projection techniques |74| . The problem of a small system weakly interacting with a heat reservoir has various aspects. 
For example, a useful model of the lattice thermal conduction is a problem of a stationary energy current through a 
crystalline lattice in contact with external heat reservoirs [TH Udi l77| . Basic to the derivation of a transport equation 
for a small system weakly interacting with a heat bath is a proper introduction of model assumptions. 
We are interested here in the problem of derivation of the kinetic equations for a certain set of average values ( 
occupation numbers, spins, etc.) which characterize the nonequilibrium state of the system. 

Let us consider the relaxation of a small subsystem weakly interacting with a thermal bath. The Hamiltonian of the 
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total system is taken in the following form: 

H = Hi + H 2 + V, (51) 

where 

H 1=Y1 S « a a a «; V = $ «/3 a l a /3> *«/3 = $ a/3 ( 52 ) 
a a,/3 

Here H\ is the Hamiltonian of the small subsystem, and a^ a and a a are the creation and annihilation second quantized 
operators of quasiparticles in the small subsystem with energies E a , V is the operator of the interaction between 
the small subsystem and the thermal bath, and H 2 is the Hamiltonian of the thermal bath which we do not write 
explicitly. The quantities $ a /3 are the operators acting on the thermal bath variables. 

We are interested in the kinetic stage of the nonequilibrium process in the system weakly coupled to a thermal bath. 
Therefore, we assume that the state of this system is determined completely by the set of averages < P a p >—< alflp > 
and the state of the thermal bath by < H 2 >, where < . . . > denotes the statistical average with the nonequilibrium 
statistical operator, which will be defined below. 

In order to pursue our discussion, we will use the whole development in section [H] We take the quasi-equilibrium 
statistical operator p q in the form 

p q (t)=exp(-S(t,0)), S{t,Q) = n(t)+Y J P a pF aP (t) + pH 2 (53) 

a/3 

n = In Tr exp(- £ P a0 F a(i (t) - f3H 2 ) 

a/3 

Here F a g(t) are the thermodynamic parameters conjugated with P a p, and /3 is the reciprocal temperature of the 
thermal bath. All the operators are considered in the Heisenberg representation. The nonequilibrium statistical 
operator has the form 



p(t)=exp(-S(t,0)), (54) 



S(t,0)=e f dhe etl I fi(i + ti) + ]T P a0 F aP (t) + (3H 2 ) 

J -°° V a/3 / 



The parameters F a p(t) are determined from the condition < P a p >=< P a p > q . 

In the derivation of the kinetic equations we use the perturbation theory in a "weakness of interaction" and assume 
that the equality < $ Q /3 > q = holds, while other terms can be added to the renormalized energy of the subsystem. 
The nonequilibrium statistical operator can be rewritten as 

p{t) = Q- 1 exp(-L(t)), (55) 

r o ( 

L(t) 



: f dt ie etl ( P a pF afi {t + 0H 2 (ti) j 
J -°° \a/J / 



Integrating in Ea. (|55|l by parts, we obtain 

L(t) = J2 P ^F a0 {t) + PH 2 (56) 



\a(3 a/3 



a/3 
dtl 



l3H 2 (h) 



For further considerations it is convenient to rewrite p q as 

Pq = P1P2 = Qg 1 exp(-L [t)), (57) 

where 

p x = Ql 1 exp P a pF aP {t) ; Q x = Trexp P a pF aP {t) (58) 

\ a/3 / \ a/3 / 
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P2 = Q^e-f^ ■ Q 2 = Tr exp(-/3P 2 ) 
Q q = Q1Q2; L = ^P a/3 F Q/3 (<) + (3H 2 



(59) 
(60) 



a/3 



We now turn to the derivation of the kinetic equations. The starting point is the kinetic equations in the following 
implicit form: 



d<P f > = L < [p ap H ] >= - E a ) < P afJ > +4 < [ p ccp, V] > 

at in in in 



(61) 



We restrict ourselves to the second-order in powers of V in calculating the r.h.s. of l|tjl|) . To this end, we must obtain 
p(t) in the first-order in V. We get 



dF a0 {t + h) _ i 



dF aP {t + h) 1 

dt x h y ~ IJ " u '" vv ' ' -i/ ^ 9<P^> £7i 



< [P^(ti),v(ti)] >-- 



~(Ep - E a )F aP (t + h) -J2 



tun 



(< $^ 7 P M7 > - < $ 7M P 7 „ >) 



(62) 



Restricting ourselves to the linear terms in Eq. l|62|) . we obtain 

dF^jt + h) i 1 dF af3 (t + h) 
— ~ -(E,3 - E a )F aP (t + h)--}^ 



(< >g< > - < $ W > g < Ry„ >) 



-(P /3 -P Q )P Q/3 (t + t 1 ) (63) 



The quantities, P a p(ti) and H 2 (ti) in the first-order in interaction have the form 

P Q/3 (il) = " E a )P a p(h) + ^[Pa/3, Vfa)] 

ff 2 (tl) = i[# 2 (*l)>^(*l)] 

Here and below all the operators are taken in the interaction representation. Using Eas. (|63[) and (|64|) we find 



L(t) = L 



dt x e E 



a(3 



(64) 



(65) 



It can be verified that the expression Y] a a Pap(ti)F a p(t + ti) + (3H 2 is independent of t\ in the zero-order in interaction 
and consequently is equal to Lq. Then for p(t) in the linear approximation in interaction V we have 



p{t) =p g - -P„ / dhe £tl / d\e XL °[L , 7(ti)]e- Ai ° 
n J -00 Jo 

By integrating in Ea. 1)66(1 over A and using the relation 



(66) 



e A ^[L ,U(<i)]e- 



-XL 



" ~ -rre XLo V(ti)e- XLo 
d\ 



we get 



p(t) = Pq --J dt ie £tl [V(h),p q } 



Finally, with the aid of Eq. (|68() we obtain the kinetic equations for < P a p > in the form 
d<P a p> 1 1 '■" 



dt 



4(P/3 ~E a ) < P aP > - i / dt ie £tl < [[P a0 ,V],V(h)] > q 

in n 7_oo 



(67) 



(68) 



(69) 
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The last term of the right-hand side of Eq.JBSJ) can be called the generalized "collision integral" . Thus, we can see that 
the collision term for the system weakly coupled to the thermal bath has a convenient form of the double commutator 
as for the generalized kinetic equations (|37|l for the system with small interaction. It should be emphasized that the 
assumption about the model form of the Hamiltonian (|51|l is nonessential for the above derivation. We can start again 
with the Hamiltonian (|51() in which we shall not specify the explicit form of Hi and V. We assume that the state 
of the noncquilibrium system is characterized completely by some set of average values < Pk > and the state of the 
thermal bath by < H2 >■ We confine ourselves to such systems for which [Hi, Pk] = c kiPi- Then we assume that 
< V > q — 0, where < . . . > q denotes the statistical average with the quasi-equilibrium statistical operator of the form 



Pq = Q- 1 exp ^- £ PkF k (t) - pH^j (70) 

and Fk(t) are the parameters conjugated with < Pk >■ Following the method used above in the derivation of equation 
(|69[) . we can obtain the generalized kinetic equations for < Pk > with an accuracy up to terms which are quadratic 
in interaction 



d<P k > 



dt 



i 1 r° 

M 1 ll J — OO 



(71) 



Hence l|69[) is fulfilled for the general form of the Hamiltonian of a small system weakly coupled to a thermal bath. 



IV. SYSTEM IN THERMAL BATH: BALANCE AND MASTER EQUATIONS 

In section UTTl we have obtained the kinetic equations for < P a p > in the general form. Our next task is to write 
down equations (|69(l in an explicit form. To do this, we note that the perturbation operator can be represented as 
y^i) = <£a/3(*i) a a a /3: where 

Kpiti) = U 2 {t 1 )^ a ^{t 1 )exp{^(E a ~ Ep)t x ); U 2 (h) = exp(?-|^) (72) 
Now we calculate the double commutator in the right-hand side of Eq. I|69|) 

< [[P a0 ,V],V(h)} > q = (73) 

^{< ®0n<f>»u(tl) >q< Potv > + < 9W(*l)$/ia >1< P "0 > ~(< ® luxfovfa) >q + < <l>»a(tl)<$>pv >q) < Pp.v >} 
//;/ 

where we restricted ourselves to the linear terms in the mean density of quasiparticles. Let us now remind that 
the correlation functions < AB(t) > and < A(t)B > can be expressed via their spectral intensities. Indeed, an 
effective way of viewing quasiparticles, quite general and consistent, is via the Green functions scheme of many-body 
theory [6^, [73. It is known H?),!!^] that the correlation functions and Green functions are completely determined by 
the spectral weight function ( or spectral intensity) J (to). 

1 / ,+oc 

F AB (t - t') =< A{t)B(t') >= — / dwexp[iw(i - t')]J AB (Lu) (74) 

1 r + °° 

FbaH' - t) =< B{t')A{t) >= — / dujexp[iu(t' - t)]J BA {u) (75) 
Here the Fourier transforms Jab{u) and Jba(^) are of the form 

Jba(u) = (76) 
Q-^^M-PEn^lB^MlA^Si ^ ~ Em -co) 

m.n 

Jab(-u) = exp(/3?iuj)JBA(u) (77) 

Expressions l|7t)|l and (|77|l are spectral representations of the corresponding time correlation functions. The quantities 
Jab and Jba are spectral densities or spectral weight functions. 
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It is convenient to define 

f +00 



F BA (0) =< B(t)A(t) >=■£-[ dtoJ(u) (78) 

2 "" J -DO 



1 



F AB (0) =< A{t)B{t) >= — / dwexp(/37iu;)J(w) (79) 

The correlation functions < ^^^(ii) > q and < </>y /J (ii) < I>^a > 9 are connected with their spectral intensities in the 
following way: 



<®uu4><y5(t)> q =— dujJ 7 s.n^(uj)exp[-i(Lj 2— )i] (80) 

/_„„ ft 



oc 
+00 



< ^ w (t)$ 7 a > g = — / rfwJ 7(5 . All/ (w)exp[i(w + — -)t] (81) 
Substituting Eas. l|80[) and (|81() into Eas. l|69[) and 173() and taking into account the notation 



ih 

"+00 



1 r 

— / *ie £tl (< $ Ma <^„(tl) >q + < ^(h^/iv >q) 



1 /"+ 00 / 1 1 \ 

— / duJp vlia (u) — — — — = K a0IMU (83) 

2tt J_ 00 m ' m \fiuj - E [3 + E u + ie hu - E a - - is J 



one can rewrite the kinetic equations for < P a p > as 
d< P a/3 > 1 



dt 



If one confines himself to the diagonal averages < P aa > only, the last equation may be transformed to give 

d < P aa > 



dt 



= K <x<x,w <Pw>- (K aa + K f aa ) < P aa > (85) 



Kaa,i/u — — -nJai/,i/a{ Z ~) — Wp—Hjt (86) 

n n 

K aa + Kl a = ^Yl J™.«A Efl h Ea ) = W a ^ (87) 

Here Wp^ a and W a — >p are the transition probabilities expressed in the spectral intensity terms. Using the properties 
of the spectral intensities [68| , it is possible to verify that the transition probabilities satisfy the relation of the detailed 
balance 

W ^ a = exp(-/3E a ) 

Finally, we have 

d < P aa > 



dt 



W v ^ a <P VV >-Y, W <*^» < P aa > (89) 



This equation has the usual form of the Pauli master equation. 

According to Ref. |79j |. "the master equation is an ordinary differential equation, describing the reduced evolution 
of the system, obtained from the full Heisenberg evolution by taking the partial expectation with respect to the 
vacuum state of the reservoirs degr ees of freedom" . The rigorous mathematical derivation of the generalized master 
equation [7^, IU |3 [III HH HH 113 is rather a complicated mathematical problem. 
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V. A DYNAMICAL SYSTEM IN A THERMAL BATH 

The problem about the appearance of a stochastic process in a dynamical system which is submitted to the 
influence of a " large" system was considered by Bogoliubov [371 H3 . For a classical system this question was studied 
on the basis of the Liouville equation for the probability distribution in the phase space and for quantum mechanical 
systems on the basis of an analogous equation for the von Neumann statistical operator. In the mentioned papers 
a mathematical method was elaborated which permitted obtaining, in the first approximation, the Fokker-Planck 
equations. Since then a lot of papers were devoted to studying this problem from various points of view ( e.g. 
Refs. |H |H |H ilil |H GjjJ). Lebowitz and Rubin H3 studied the motion of a Brownian particle in a fluid ( 
as well as the motion of a Brownian particle in a crystal ) from a dynamical point of view. They derived a formal 
structure of the collision term similar to the structure of the usual linear transport equation. Kassner |88| used a 
new type of projection operator and derived homogeneous equations of motion for the reduced density operator of a 
system coupled to a bath. It was shown that in order to consistently describe damping within quantum mechanics, 
one must couple the open system of interest to a heat reservoir. The problem of the inclusion of dissi pat ive forces in 
quantum mechanics is of great interest. There are various approaches to this complicated problem |79l \83L l9ll 192 . l93| . 
Tanimura and Kubo [9l| considered a test system coupled to a bath system with linear interactions and derived a set 
of hierarchical equations for the evolution of their reduced density operator. Breuer and Petruccione developed 
a formulation of quantum statistical ensembles in terms of probability distributions on a projective Hilbert state. 
They derived a Liouville master equation for the reduced probability distribution of an open quantum system. It 
was shown that the time-dependent wave function of an open quantum system represented a well-defined stochastic 
process which is generated by the nonlinear Schrodinger equation 

^ = -iGty) (90) 

with the nonlinear and non-Hermitian operator G(ip). The inclusion of dissipative forces in quantum mechanics 
through the use of non-Hermitian Hamiltonians is of great interest in the theory of interaction between heavy ions. 
It is clear that if the Hamiltonian has a non-Hermitian part Ha the Heisenberg equation of motion will be modified 
by additional terms. However, care must be taken in defining the probability density operator when the Hamiltonian 
is non-Hermitian. Also, the state described by the wave function ip is not then an energy eigenstate because of the 
energy dissipation. As it was formulated by Accardi and Lu [79j |. "the quantum Langevin equation is a quantum 
stochastic differential equation driven by some quantum noise ( creation, annihilation, number noises)." The necessity 
of considering such processes arises in the description of various quantum phenomena ( e.g., radiation damping, etc.), 
since quantum systems experience dissipation and fluctuations through interaction with a reservoir |94l l95| . The con- 
cept of "quantum noise" was proposed by Senitzky |93j to derive a quantum dissipation mechanism. Originally, the 
time evolution of quantum systems with the dissipation and fluctuations was described by adding a dissipative term 
to the quantum equation of motion. However, as was noted by Senitzky :93], this procedure leads to the nonunitary 
time evolution. He proposed to derive the quantum dissipation mechanism by introducing quantum noise, i.e., a 
quantum field interacting with the dynamical system (in his case an oscillator). For an appropriately chosen form of 
the interaction, energy will flow away from the oscillator to the quantum noise field (thermal bath or reservoir). 
In this section, we consider the behavior of a small dynamic system interacting with a thermal bath, i.e., with a 
system that has effectively an infinite number of degrees of freedom, in the approach of the nonequilibrium statistical 
operator, on the basis of the equations derived in section ITTT1 The equations derived below can help in the under- 
standing of the origin of irreversible behavior in quantum phenomena. 

We assume that the dynamic system ( system of particles) is far from equilibrium with the thermal bath and cannot, 
in general, be characterized by a temperature. As a result of the interaction with the thermal bath, such a system 
acquires some statistical characteristics but remains essentially a mechanical system. Our aim is to obtain an equation 
of evolution ( equations of motion ) for the relevant variables which are characteristic of the system under consid- 
eration. The basic idea is to eliminate effectively the thermal bath variables (c.f. Ref. 0, |9f|, |9(j ). The influence 
of the thermal bath is manifested then as an effect of friction of the particle in a medium. The presence of friction 
leads to dissipation and, thus, to irreversible processes. In this respect, our philosophy coincides precisely with the 
Lax statement 94] "that the reservoir can be completely eliminated provided that the frequency shifts and dissipation 
induced by the reservoir are incorporated into the mean equations of motion, and provided that a suitable operator 
noise source with the correct moments are added" . 

Let us consider the behavior of a small subsystem with Hamiltonian H\ interacting with a thermal bath with Hamil- 
tonian i?2- The total Hamiltonian has the form l|51f). As operators P m determining the nonequilibrium state of the 
small subsystem, we take a^,a Q , and n a = a) a a a . Note that the choice of only the operators n a and Hi would lead 
to kinetic equations l|71|l for the system in the thermal bath derived above. 

The quasi-equilibrium statistical operator JSJ is determined from the extremum of the information entropy Q subject 
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to the additional conditions that the quantities 

Tr(pa a ) =< a a >, Tr(pa ] a ) =< >, Tr(pn a ) =< n a > (91) 
remain constant during the variation and the normalization Tr{p) = 1 is preserved. The operator p q has the form 

Pq = eiq ? ^-J2(U(t> a + fUt)ai+F a {t)n a )-f3H^ = exp(S(t,0)) (92) 

fi = bxTr exp ^- £(/a(*)a Q + /+(i)a+ + F Q (i)n Q ) - /3if 2 l 

Here, f a , /t and i 7 ^ are Lagrangian multipliers determined by the conditions l|91|l . They are the parameters conjugate 
to < a a > q , < aj, > q and < n a > q : 

ss , , . 55 , . 

< Q a > g = - 77 , , , < n a > = - — — = f a {t), — — = F Q t 93) 

<5/o (*) oF Q (t) S < a a > q 6<n a > q 

In what follows, it is convenient to write the quasi-equilibrium statistical operator (|92() in the form 

Pq=PlP2, (94) 

where 



pi = exp - 5^(/ a (t)oa + /<£(*)4 + F Q (<)n Q )^J 
fil =lnTrcxp ^- ^(/ Q (t)a a + /j^a* + F a (t)n a )^ 



(95) 



p 2 = exp (fi 2 -0H 2 ), n 2 = In Tr exp (-/3if 2 ) (96) 
The nonequilibrium statistical operator p will have the form l|23l) . Note, that the following conditions are satisfied: 

< a a > q =< a a >, < > q =< a* >, < n a > q =< n a > (97) 

We shall take, as our starting point, the equations of motion for the operators averaged with the nonequilibrium 
statistical operator (|23|l 

iK d < °° > =< [o a ,ffi] > + < [a a ,V] >, (98) 
d <c n > 

i?i -2— =< [n a , Hi] > + <[n a ,V]> (99) 

at 

The equation for < a' a > can be obtained by taking the conjugate of (|98|) . Restricting ourselves to the second order 
in the interaction V, we obtain, by analogy with l|71|l . the following equations: 

d<a a > „ ^ _ , 1 f° „ rt 



in -2— = £ Q < a Q > +— / di ie £tl < [[a a ,V], Vih)} > q (100) 

ai in ' 



—00 




ih d< 2 a> = i / * ie £tl < [K, v], >, (101) 
df »nj_ oc 

Here V(ii) denotes the interaction representation of the operator V. Expanding the double commutator in Ea. HlOOfl . 
we obtain 



= E a < a a > +i J dt ie Etl f < $«/3?W(*l) >g< 



ih ~~dT~ = E a <a a > / dtie etl ( ^ < $ a p4>iw{ti) >q< apa\a v > q - < ^(*i)$ a/ 3 > g < a\a u a/3 > 5 



(102) 
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where M „(ii) = ®nv(h) e-np(j-(E fl - E v )t\). We transform Ea. (|102[l to 
d <C (2 ^> 1 /* 

3T— = < a Q > +TT- ^ / eftie £tl < ^c^pih) > q < ap > + (103) 

(It «ft ^ ./_«, 



1 f 

" a J — OO 



We assume that the terms of higher order than linear can be dropped in Eq. l|103[) ( below, we shall formulate the 
conditions when this is possible). Then we get 

ih - < "° > = E a < a a > +i Y [ dt ie etl < Qanfapih) > q < ap > (104) 
at m ^ ./_«, 

The form of the linear equation (|104f> is the same for Bose and Fermi statistics. 

Using the spectral representations, Ea. (|80|l and Ea. (|81|l . it is possible to rewrite Ea. (|104|) by analogy with Ea. (|84fl as 

ih - < a ° > = E a < a a > + Y K a() <a f3 > (105) 
dt 

where K a p is defined in l|82l) . Thus, we have obtained the equation of motion for the average < a a >. It is clear that 
this equation describes approximately the evolution of the state of the dynamic system interacting with the thermal 
bath. The last term in the right-hand side of this equation leads to the shift of energy E a and to the damping due 
to the interaction with the thermal bath ( or medium ). In a certain sense, it is possible to say that Ea. (|105|l is an 
analog or the generalization of the Schrodinger equation. 

Let us now show how, in the case of Bose statistics, we can take into account the nonlinear terms which lead to a 
coupled system of equations for < a a > and < n a >. Let us consider the quantity < aj^auap > q . After the canonical 
transformation 

a a = b a + < a a >, a a = b a + < a a > 
the operator pi in Ea. l|95|l can be written in the form 



pi=Qr 1 expfni-^(F a (t)6t6 a )j , <aa>= -A (106) 

Note that Qi in (|106fl is not, in general, equal to Qi in (|95|l . Using the Wick-De Dominicis theorem (6^1 f° r the 
operators b a , b a and returning to the original operators a a ,a a , we obtain 

< a^a v ap > q ~ (< > -| < a M > | 2 ) < a„ > d^p + (< > -| < a M > | 2 ) < ap > S^p (107) 

Using H107J1 . we can rewrite Ea. (|95|) in the form 

d < a a > 



1 f° 

E a <a a > +— Y / dhe Etl < ^a^pih) > q < a fi > + 
Pp. J -°° 

1 f° 

/ dtie£tl (< <?W(*l)] >q + < [®cxP, >q) (< n M > +1 < a M > I 2 ) < a > 

l ' 1 . a J —OO 



(108) 



i. . 

pP 



Now consider Ea. (|lt)l(l . Expand the double commutator and, in the same way as the threefold terms were neglected 
in the derivation of Ea. H104f) . ignore the fourfold terms in l|101|) . We obtain then 

d < "° > = Y W 0^{< np>+\<ap> | 2 ) - W a ^p(< n a >+\<a a > | 2 ) + (109) 

Y Ka P < a a >< a > Y K aP < a a >< dp > + ^ K aa ^ v < a* >< a v > 

p p p.v 
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Thus, in the general case Eas. (|100|) and (|101|) form a coupled system of nonlinear equations of Schrodinger and kinetic 
types. The nonlinear equation (|102fl of Schrodinger type is an auxiliary equation and, in conjunction with the equation 
of kinetic type (|109|l . determines the parameters of the nonequilibrium statistical operator since in the case of Bose 
statistics 

<a >= <n >= ( e F a (t) _ -n-i , i/"i 2 n 101 



Therefore, the linear Schrodinger equation is a fairly good approximation if 



(< n a > +| < a a > | 2 ) = (e F °W - l)" 1 



< 1 



The last condition corresponds essentially to < b^ba 1. 

In the case of Fermi statistics the situation is more complicated 1971. There is well-known isomorphism between 
bilinear products of fermion operators and the Pauli spin matrices |98|. In quantum field theory the sources linear 
in the Fermi operators are introduced by means of classical spinor fields that anticommute with one another and 
with the original field. The Fermion number processes in the time evolution of a certain quantum Hamiltonian model 
were investigated in Ref. [9fj. It was shown that the time evolution tended to the solution of a quantum stochastic 
differential equation driven by the Fermion number processes. We shall not consider here this complicated case. 
In order to interpret the physical meaning of the derived equations, an example will be given here. Let us consider 
briefly a system of electrons in a lattice described by the Hamiltonian 

H = H 1+ H 2 + V = ^2 <k)ala ka + ]T Hu q b\b q + ± £ A(£ - k^a^b^ + b^_.J, (111) 

kcr q ki,k2(7 

where fvuj q is the phonon energy, a klr ,akcr and b q ,b q are the operators of creation and annihilation of electrons and 
phonons, respectively; e(k) is the energy of electrons and A(q) determines the electron-phonon coupling. Equation 
(|105fl for < a ka > can be represented in the form 

* h d< 2 a> = (e(fc) + AE ( k)) < ak ° > -Y T{k) < aka >J (112) 



where 



AE{k) = PJ2 ki)\ 2 f " 

fci V' 



< N k _ kl > +1 < N k - kl > 



(113) 



e(k) - e(fci) - e(k) - e(fci) + f^%_ kl 

r ( fc ) = \ A ^- k ^ 2 ((< N *-* > + 1 ) S « k ) - < k ^) ~ < N k-k! > S(e(k) - efa) + ^_ fcl ))(114) 



are the energy shift of an electron and the electron damping, respectively. Here < N q >= (^e l3huJq — 1) 1 , the 
distribution functions of the phonons. Expressions 1)1131) and l|114|) are the same as those obtained by the Green 
functions method 100] if one sets < a ka a ka ><C 1 in the latter. 



VI. SCHRODINGER-TYPE EQUATION WITH DAMPING FOR A DYNAMICAL SYSTEM IN A 

THERMAL BATH 

In the previous section we obtained an equation for mean values of the amplitudes in the form l|105|) ■ It is of interest 
to analyze and track more closely the analogy with the Schrodinger equation in the coordinate form. To do this, by 
convention we define the " wave function" 

a 

where {Xa(r)} is a complete orthonormalized system of single-particle functions of the operator ^— ^-V 2 +v{fu, 
where v(r) is the potential energy, and 

h 2 \ 

— V 2 + v(r) Xa (r) = E aXa {r) (116) 
2m / 
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Thus, in a certain sense, the quantity ip(f) may plays the role of the wave function of a particle in the medium. Now, 
using l|115|) . we transform Eq. l|lU5fl to 

„. fc W) = f_^_y2 +v ^\ ^ + f K(r,r')ip(r')d? (117) 



dt \ 2m 

The kernel K(r, r') of the integral equation l|117[) has the form 

K(f, ?) = J2 K apXa (r)xW') =^J2 [ dt ^ £tl < >, Xa(f)x},{?) (H8) 

Equation l|117|) can be called a Schrodinger-type equation with damping for a dynamical system in a thermal bath. 
It is in teresting to note that similar Schrodinger equations with a nonlocal interaction are used in the scattering 
theory jlOlL Il02fl to describe interaction with many scattering centers. 

To demonstrate the capabilities of equation (|117|l . it is convenient to introduce the operator of translation exip(iqp/h), 
where q = r' — f;p = —iKS/ r . Then Eci. (|117|) can be rewritten in the form 

m dMp_ = f±_ V 2 + v(rl \ m + J2D(r,mr) (119) 
\ v 

where 

D(r,p) = J d 3 qK{r,f+q)e^ (120) 

It is reasonable to assume that the wave function i/i(r) varies little over the correlation length characteristic of the 
kernel K(f,r'). Then, expanding cx^{iqp/K) in a series, we obtain the following equation in the zeroth order: 

ih dW)_ = f_iLv 2 + v (r + ReU(r))\ ip(r) + HmU(r)ip(r) (121) 
at \ 2m J 

where 



U(r) = ReU(r) + ilmU(r) = J d 3 qK{r, f+ q) (122) 

Expression (|121fl has the form of a Sch rodinger equation with a complex potential. Equations of this form are 
well known in the scattering theory 102] in which one introduces an interaction describing absorption (Im[/(r) < 0). 
Further, expanding exjp(iqp/h) in a series up to the second order inclusively, we can represent Ea. l|117| in the following 
form [l02]: 

ih <^ = { f_tf w 2 +v ^\ +uif) _l I SqK^+^^ + l j d 3 qK(r,r + q) ]T <z" l g"V m V„}^(r) (123) 
\ m / 1 J J m,n=l 

To interpret this equation, let us introduce the function 

Aif) = ^ / d 3 qReK(r,r + q)q (124) 
ine J 

where m and e are the mass and charge of the electron and c is the velocity of light. Then A(r) can be considered, in 
a certain sense, as an analog of the complex vector potential of an electromagnetic field. It is clear that the motion 
of a particle (dynamic subsystem) through the medium imitates, to some extent, the motion of a charged particle in 
the electromagnetic field. To make this analogy even more close, let us introduce the following quantity: 

I \ Til C I 

1 = 3% ~ — / d q)q l q> (125) 



M (r) / i ■ m ihe 

It follows from l|125|) that this quantity can be int erpre ted as a tensor of the reciprocal effective masses |l03l Il04| . 
The notion of the "mass tensor" was introduced in |l03l | to describe the motion of an electron in an external field F 

dvi e v-^ d 2 E 

-r- = ^r> „ - Fi i, j ' = 1,2, 3 or x,y,z 126 

dt h 2 ^ dq l dqJ ' w ' y ' 

3 
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or m vector notation 

dv e 



— = ^grad^Fgrad^) (127) 



Thus, a field F may change the velocity v in directions other than that of F. The quantity Ti 2 ( d d ql g q j ) 1 has been 
called the "mass tensor". Now we can rewrite Ea. (|123[l in the form 

T(r) = ^ f d 3 qK(r,r + q) ^ g"Y 1 V„ l V„ (129) 

^ m,n=l 

Note that in an isotropic medium the tensor f ^ is diagonal and A(r) = 0. The introduction of ip(r) does not 

mean that the state of the small dynamical subsystem is pure. It remains mixed since it is described by the statistical 
operator l|2^ll , the evolution of the parameters f a , fl, and F a of the latter being governed by a coupled system of 
equations of Schrodinger and kinetic types. It is interesting to mention that the derivation of a Schrodinger-type 
equation with non-Hcrmitian Hamiltonian which describes the dynamic and statistical aspects of the motion was 
declared by Korringa jg^. However, his Eq.(29) 

dW ( i dh! \ 

t^-= ( H \t)+h'(t)) + ±=2- + ...)w'{t) (130) 



dt V 29 dt 

where W'(t) is the statistical matrix for the primed system, can hardly be considered as a Schrodinger-type equation. 
This special form of the equation for the time-dependent statistical matrix can be considered as a modified Bloch 
equation. 

Hence we were able to apply the NSO approach given above to dynamics. We have shown in this section that for some 
class of dynamic systems it was possible, with the NSO approach, to go from a Hamiltonian description of dynamics 
to a description in terms of processes which incorporates the dissipativity. However, a careful examination is required 
in order to see under what conditions the Schrodinger-type equation with damping can really be used. 



VII. CONCLUDING REMARKS 



In this paper, we have discussed the general statistical mechanics approach to the description of the transport 
processes. We have applied the method of the nonequilibrium statistical operator to study the generalized kinetic 
and evolution equations. We analyzed and derived in a closed form the kinetic equations and applied them to some 
typical problems. 

In writing the paper we have essentially confined ourselves to a discussion of those features of the theory which deal 
with general structural properties rather than with specific physical applications. The method offers several advantages 
over the standard technique of the calculation of transport coefficients. The derived generalized kinetic equations for 
a system weakly coupled to a thermal bath are analogous to those derived in for the system of weakly interacting 
particles. Moreover, the capability of the generalized kinetic equations was demonstrated and further discussed by 
considering a few representative examples, i.e., the kinetic equations for magnons and phonons, and the energy shift 
and damping of particle (electron) due to the friction with media (phonons). There are many other applications of 
the formalism developed in this article, for example, longitudinal nuclear spin relaxation and spin diffusion. However, 
we have not considered other contributions here. These questions deserve a separate consideration. 
An example of a small system being initially far from equilibrium has been considered. We have reformulated the 
theory of the time evolution of a small dynamic system weakly coupled to a thermal bath and shown that a Schrodinger- 
type equation emerges from this theory as a particular case. Clearly then, the nonequilibrium statistical operator 
approach is a convenient and workable tool for the derivation of relaxation equations and formulae for evolution and 
kinetic equations. 

In our above treatment we have avoided a number of important questions such as the rigorous proof of the existence 
and uniqueness of the quasi-equilibrium state, the validity of the time-smoothing procedure, etc. These questions, as 
well as the application of the derived equations to other important problems of transport in solids such as the nuclear 
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spin relaxation and diffusion, electro- and thermal conductivity, remain to be areas for further investigation. 
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